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We  report  precise  measurements  of  the  enthalpy  of  polycrystalline 
YBa2CU307.3  at  high  temperatures  (300  K  <  T  <  700  K)  .  From  these 
data,  we  determine  the  lattice  and  electronic  contributions  to 
the  specific  heat.  We  find  that  anharmonicity  must  be  included 
when  calculating  the  lattice  specific  heat.  We  also  find  that 
the  specific  heat  of  the  charge  carriers  follows  the  temperature 
dependence  of  a  Fermi  liquid  with  a  Sommerfeld  constant  of  40±5 
mJ/mol»K^.  This  value  is  three  times  larger  than  expected  from 
band  structure  calculations. 
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Much  effort  has  been  focused  on  experimentally  determining 
the  value  of  the  coefficient  of  the  electronic  specific  heat, 

for  YBa2Cu307-4 .  An  experimental  measure  of  yg  is  needed  to  test 
the  theoretical  models  for  the  band  structure  in  YBa2Cu307.^. 

Also,  the  value  of  this  parameter  is  needed  to  determine  the 
relative  size  of  the  specific  heat  jump  at  so  that  theoretical 
models  of  the  superconducting  state  can  be  tested.  Despite  the 
importance  of  and  the  level  of  this  effort,  the  value  of  y^  must 
be  judged  to  be  uncertain  as  estimates  range  from  13  to  35 
mJ/mol»K^^-2 

Several  groups  have  attempted  to  determine  yg  from  specific 
heat  data  measured  between  the  superconducting  transition 
temperature  (T^)  and  300  However,  the  electrons  contribute 

only  1%  to  the  specific  heat  in  this  temperature  range,  with  the 
remainder  due  to  the  lattice.  Therefore,  the  lattice 
contribution  would  have  to  be  calculated  with  an  accuracy 
difficult  to  attain  in  order  to  determine  the  electronic  specific 
heat  from  normal  state  measurements  below  300  K. 

At  high  temperatures  (T  >  300  K) ,  the  lattice  specific  heat 
approaches  the  classical  limit,  while  the  electronic  contribution 
continues  to  increase  linearly.  Thus,  the  possibility  for 
accurately  modeling  the  lattice  contribution  improves,  and  the 
electronic  term  grows  into  a  larger  fraction  of  the  total.  Two 
groups*  have  measured  the  specific  heat  of  YBa2Cu307.i  at  high 
temperatures  using  differential  scanning  calorimetry.  This 
method  lacks  the  necessary  accuracy  to  determine  the  electronic 
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specific  heat.  We  have  therefore  used  a  phase-change  calorimeter 
to  measure  the  temperature  dependence  of  the  enthalpy  of 
polycrystalline  YBa2CU307.^  between  300  K  and  700  K.  The  data  are 
sufficiently  accurate  that  the  electronic  contribution  can  be 
determined  by  subtracting  the  lattice  specific  heat  (c^*^) 
calculated  from  the  phonon  density  of  states.  We  shall  show  that 
the  role  of  anharmonicity  is  significant  in  the  lattice  specific 
heat,  and  the  electronic  contribution  is  larger  than  previously 
estimated . ^ 


The  structure  of  YBa2Cu307.f  forms  two  subsystems:  CUO2  planes 


which  extend  in  the  a  and  b  directions  and/ CuO ^hains  which 

- '4-i^ 

follow  the  b-axis.^  The  oxygen  on  the  chain-sites''^ is  much  more 
volatile  than  on  the  plane-sites,  so  the  oxygen  deficiency,  6,  is 
mostly  due  to  vacancies  in  the  chains.  The  electronic  properties 
are  observed  to  be  extremely  sensitive  to  the  oxygen 
stoichiometry.  Hence,  we  have  been  careful  to  prepare  samples 
with  a  small  value  of  ft  and  to  measure  the  specific  heat  under 
such  conditions  that  oxygen  is  not  lost  during  the  course  of  the 
experiment . 

Polycrystalline  YBa2Cu307.3  was  prepared  by  solid-state 
reaction  in  air  of  stoichiometric  quantities  of  Y2O3,  Ba02,  and 
CuO  to  yield  approximately  40  g  of  product.  Reactants  were  mixed 
before  each  heating  by  grinding  for  50  minutes  in  an  agate  mortar 
and  pestle.  Two  heating  cycles  were  carriel  out  in  an  alumina 
combustion  boat:  1)  Heat  to  880*C  in  3h,  anneal  24h,  furnace 


cool;  2)  Heat  to  950*C  in  3h,  anneal  15h,  cool  in  8h  to  650'C, 
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cool  in  8h  to  450'C,  anneal  12h,  cool  in  8h  to  ambient.  The 
fully  reacted  powder  was  then  annealed  at  450“C  for  24h  in 


flowing  oxygen.  The  product  was  '^,pure  phase  according  to  X-ray 

o 

powder  diffraction,  all  observed  reflections  matching  those  of 
the  reported  powder  pattern.®  Least  squares  refinement  of  the 


powder  data  yielded  the  following  lattice  parameters:  a  = 
3.8838(8),  b  =  3.8164(4),  and  c  =  11.664(1)A.  Bromo/iodimetric 
titration^  of  the  powder  indicated  an  oxygen  content  of  6.98  per 
mole  compound  (i.e.,  6  =  0.02).  DC  magnetic  susceptibility 
measurements  showed  the  onset  of  superconductivity  at  92  K  with  a 
10  to  90%  transition  width  of  6  K.  We  observed  that  45%  of  the 
flux  was  expelled  when  the  sample  was  cooled  from  100  K  to  4.2  K 
in  a  0.002  tesla  field. 

The  method  for  measuring  specific  heat  by  means  of  phase- 
change  calorimetry  is  described  fully  in  Refs.  8.  The  sample 
(mass:  4.7g)  is  encapsulated  with  exchange  gas  in  a  Pt-10%Rh 
vial.  (We  chose  33  kPa  of  O2  at  300  K  so  that  the  oxygen  content 
of  the  sample  would  not  change  when  the  sample  was  heated.)  The 
encapsulated  sample  is  heated  to  thermal  equilibrium  in  a 
furnace,  and  then  dropped  into  a  calorimeter  immersed  in  a  high 
purity  ice  bath.  The  enthalpy  released  by  cooling  the  sample 
from  T  to  273.15  K  melts  a  proportional  volume  of  ice.  The 
relative  enthalpy  is  then  calculated  from  the  change  in  volume  of 
the  water-ice  system  employing  an  independent  electrical 
calibration.  The  measured  values  of  the  relative  enthalpy  of 
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YBajCUjOj-j  (Hp(t)  -  Hp(0))  are  plotted  vs.  temperature  ('C)  in 
Fig.  1. 

We  note  that  this  method  accurately  measures  the  enthalpy 
between  a  temperature  T  and  273.15  K.  The  sample  temperature  is 
known  to  within  10  mK  and  enthalpy  changes  are  measured  with  an 
imprecision  of  0.5  J.  The  enthalpy  of  the  empty  cell  is  measured 
in  a  second  experiment.  For  the  measurements  reported  here,  the 
addenda  contribute  half  of  the  total  enthalpy  over  the  entire 
temperature  range.  The  estimated  overall  inaccuracy  of  an 
individual  enthalpy  datum  is  0.5%.  Specific  heat  data  derived 
from  the  enthalpy  are  estimated  to  have  an  inaccuracy  not 
exceeding  1%. 

We  calculate  the  average  specific  heat  (5p  ■  AHp/t)  for  each 
datum  of  Fig.  1  and  fit  these  values  with  the  function  a  +  bt  + 
ct^,  where  t  is  the  Celsius  temperature.  The  true  specific  heat 
values  are  calculated  from  the  relation, 

Cp  •  Cp  +  t{dc~p/dt)  =  dAHp/dt  =  AHp/t  *  bt  *  2ct^  (1) 

Finally,  we  calculate  the  specific  heat  at  constant  volume  (Cy) 
by  subtracting  the  dilatation  correction,  P^BVT  (P  is  the  volume 
thermal  expansion  coefficient  (36  x  10  ®  K'^)®,  B  is  the  bulk 
modulus  (10^^  Pa)^°,  V  is  the  molar  volume  (1.04  x  10'‘  mVniol)^, 
and  T  is  the  absolute  temperature.  (The  temperature  dependence 
of  P  and  V  has  been  shown  to  be  small  between  300  K  and  700  K,“ 
but  the  temperature  dependence  of  the  bulk  modulus,  which  we 
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assume  to  be  constant,  has  not  been  measured  above  300  K.)  The 
resulting  values  for  the  specific  heat  at  constant  volume  are 
plotted  in  Fig.  2.  Note  that  the  specific  heat  data  measured 
below  300  (smaller  circles)  are  in  good  agreement  with  the 
higher  temperature  data. 

The  coefficients  of  Eg.  1  have  been  tested  for  uniqueness  by 
fitting  H  vs.  t  to  a  polynomial,  at  +  bt^  +  ct^.  We  find 
slightly  different  values  for  a,  b,  and  c,  because  more  weight  is 
given  to  the  higher  temperature  points  in  the  fits  of  H  vs.  t 
than  in  the  fits  of  H/t  vs.  t.  Nevertheless,  the  difference 
between  the  c^  vs.  T  curves  calculated  from  the  two  sets  of  fit 
parameters  is  slight. 

We  calculate  the  lattice  contribution  to  c^  from  the  phonon 
density  of  states,  F((i)),  determined  from  the  measured  neutron 
density  of  states,  (o)  .  ^  We  calculate  the  harmonic 

contribution  to  the  heat  capacity  by  directly  integrating  the 
contribution  of  each  mode  in  the  phonon  spectrum  according  to  the 


formula , 


c 


j 

V 


j"da)  F  ( 0) )  o)  f  ^  ( CO  ,  T) 

,0  t 


(2) 


where  fB((o,T)  is  the  Bose-Einstein  distribution.  When  this 
lattice  component  is  subtracted,  we  are  left  with  the  electronic 
contribution  to  the  specific  heat  seen  in  the  inset  of  Fig.  3. 

If  the  charge  carriers  pair  to  form  bipolarons  in  the  normal 
state,”  we  would  expect  the  electronic  specific  heat  to  be 
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similar  to  that  of  a  Boson  gas.  That  is,  the  specific  heat 
should  peak  at  the  condensation  temperature  (T^.  =  92  K  in  the 
case  of  YBa2Cu307.^)  and  then  decrease  monotonically  to  the 
equipartition  value.  Our  data  are  therefore  inconsistent  with 
this  model.  As  shown  in  the  inset  of  Fig.  3,  the  temperature 
dependence  of  the  electronic  contribution  is  characteristic  of 
the  specific  heat  of  a  normal  Fermi  liquid,  except  that  a  line 
fitted  to  the  data  does  not  pass  through  the  origin.  The 
negative  intercept  implies  that  the  contribution  of  Eq.  2 
overestimates  the  lattice  specific  heat,  probably  due  to  the 
effect  of  lattice  anharmonicity .  Theoretical  calculations,^^ 
Mossbauer  effect  measurements,^®  and  the  anomalous  pressure 
dependence  of  the  compressibility^®  all  show  evidence  of 
anharmonicity  in  the  zone-center  phonons.  This  suggests  that  the 
harmonic  approximation  is  inadequate. 

The  energy  of  an  atom  in  a  periodic  potential  well  can  be 
written  E  =  pVsm  +  aq"  (n  is  an  even  integer  greater  than  or 
equal  to  zero,  n  =  2  for  a  harmonic  potential)  .  Hui  and  Alien®® 
have  shown  that  a  quartic  anharmonicity  reduces  the  high- 
temperature  specific  heat  relative  to  the  value  expected  from  a 
harmonic  potential.  In  fact,  the  equipartition  value  for  the 
specific  heat  per  mode  is  kB(l/2  +  l/n)  for  an  atom  in  a 
potential,  aq".  Clearly,  the  specific  heat  for  purely  harmonic 
modes  will  reach  the  Dulong-Petit  limit  (ks  per  mode).  However, 
the  Dulong-  Petit  value  overestimates  the  contribution  to  the 
specific  heat  from  terms  in  the  interatomic  potential  with  n  >  2. 
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Since  the  exact  nature  of  the  anharmonicity  in  YBa2Cu307,^  is 
unknown  at  this  time,  we  propose  a  simple  model  that  accounts  for 
the  anharmonicity  in  a  way  that  produces  a  physically  reasonable 
electronic  specific  heat:  the  interatomic  potential  for  a 
fraction,  n,,  of  the  phonons  with  energies  less  than  20  meV  is 
quart ic.  Thirty  percent  of  the  phonons  lie  below 
20  meV  and  these  will  have  reached  the  equipartition  limit 
(0.3x3Rx(l/2  +l/n) ,  R  is  the  universal  gas  constant)  by  the 
temperature,  300  K.  The  total  specific  heat  for  temperatures 
greater  than  300  K  can  then  be  written; 


=  0.9i?(l  - 


+  0.675Rn.  -t- 


dT 


/ 


dti)F(a>)  wfu  (g)  ,  r) 


Vo. 02 A 


(3) 


Clearly,  as  n,  is  increased,  the  lattice  specific  heat  will  be 
reduced. 

By  varying  n,  until  the  electronic  specific  heat 
extrapolates  to  the  origin,  we  determine  the  electronic  specific 
heat  plotted  in  Fig.  3.  If  the  upper  and  lower  lines  bracket  the 
electronic  specific  heat,  we  find  that  n,  is  0.25±0.08  and  the 
slope,  Ysf  is  40±5  mJ/mol-K^.  This  is  higher  than  other  reported 
values^,  the  highest  of  which  are  in  the  range,  20  to  35 
iiiJ/mol»K^.  These  estimates  for  Yb  determined  from  specific 
heat  data  measured  below  300  K;  a  temperature  range  in  which  the 
lattice  specific  heat  is  difficult  to  model. 

The  measured  Sommerfeld  constant,  40  mJ/mol*K^,  differs  from 
the  value  of  13  mJ/mol«K^  obtained  from  band  structure 
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calculations.^^  This  discrepancy  could  be  due  to  the  complicated 
band  structure  in  YBa2Cu307.^.  The  Fermi  energy  lies  near  the  foot 
of  a  peak  in  the  calculated  electronic  density  of  states,  D(Cp).*^ 
Hence,  sample  stoichiometry,  defects,  or  inaccuracies  in  the 
calculation  could  shift  the  electronic  structure  enough  to  raise 
D(ep)  to  a  value  more  consistent  with  our  measurements. 

The  implications  of  a  large  value  for  Yb  best  understood 

in  terms  of  a  two  band  model.  The  resistivity  of  untwinned 
single  crystals^^  shows  the  existence  of  two  conducting  bands: 
one  attributed  to  the  chains  and  the  other  to  the  planes.  These 
same  data  show  that  the  major  fraction  of  the  current  is  carried 
by  the  chains.  Other  studies^*  have  shown  that  as  5  increases 
from  0  to  0.2,  T^  remains  nearly  constant  and  the  valence  bond 
sum  of  the  plane  copper  atoms  changes  very  little.  Thus,  when  a 
small  amount  of  oxygen  (6  s  0.2)  is  removed,  the  electronic 
properties  of  the  chains  are  more  strongly  affected  than  those  of 
the  planes.  Evidence  that  the  chain  band  contributes  greatly  to 
D(ep)  is  seen  in  specific  heat  measurements  of  YBa2CU307-j  with  6  = 
0.02  and  6  =  0.2.'^^  The  difference  between  the  specific  heat  of 
these  two  materials  is  linear  in  temperature  with  a  slope,  32 
mJ/mol*K^.  By  identifying  slope  of  this  difference  with  the 
electronic  specific  heat  of  the  chain  band,  we  utilize  our 
maximum  value  of  Yb  (i-e.  Ych  +  Ypi)  to  estimate  that  Ypi  ^  13 
mJ/mol*K^.  Thus  the  plane-carrier  effective  mass  (in  two 
dimensions,  m*pi  =  ( 3/ti )  (Yi/ke)  ^dYpi,  where  d  is  the  interplanar 
distance)^®  is  less  than  9m,  (m,  is  the  f ree-electron  mass). 
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The  ratio  of  the  Soromerfeld  constants  for  the  plane  aid 
chain  bands  affects  the  interpretation  of  the  temperature 
dependence  of  the  penetration  depth  and  the  infrared  (ir) 
reflectivity.  If  each  band  is  characterized  by  the  ratio  of  the 
carrier  density  to  the  effective  mass,  n/m" ,  then  the  relative 
size  of  this  ratio  is  defined  to  be  Rq  =  (npi/rr*pi)  /  (nch/m%h)  •  We 
use  the  expressions  for  the  effective  mass  of  a  two  dimensional 
band  (see  previous  paragraph)  and  of  a  one  dimensional  bana  (m^j, 

=  12  (lidb/kg)  ^n^hYch)  to  show  that  Rq  =  47idb^npiYch/Ypi  is  the 
interplanar  spacing  (11. eA)^,  b  is  the  distance  between  chains 
(l.SsA)^,  npi  is  the  density  of  carriers  in  vhe  planes  (2.86x10^^ 
cro’^)^®).  Since  the  planes  remain  stable  for  oxygen  deficiencies 
less  than  6  =  0.2,  we  contend  that  Rq  is  affected  mostly  by  the 
structural  integrity  of  the  chains.  The  specific  heat  data 
indicate  that  Rq  is  greater  than  or  equal  to  15,  for  fully 
oxygenated  samples.  This  compares  to  the  value  of  3.5  determined 
from  penetration  depth  measurements^^  in  thin  films  and  0.5  from 
the  ir  reflectivity  of  untwinned  single  crystals.^®  Noting  that 
the  penetration  depth  data  of  Anlage  et  al.^^-  are  strongly 
dependent  on  the  number  of  defects  in  the  film,  we  expect  that 
films  with  fewer  chain  defects  would  exhibit  a  higher  vaj^e  of 
Rq.  Similarly,  the  value  of  Rq  determined  from  ir  reflectivity 
data  should  increase  in  samples  with  higher  oxygen  stoichiometry. 
In  general,  we  speculate  that  the  strong  dependence  of  the 
measured  electronic  properties  of  YBaiCUjO;.^  (i.e.  specific  heat 


If 


jump  at  T^,  penetration  depth,  and  resistivity)  on  oxygen  content 
is  due  to  the  large  contribution  of  the  chain  band  to  D(ef)  . 

The  high  temperature  specific  heat  reveals  much  about  the 
electronic  and  lattice  dynamical  properties  of  YBa2Cu307.^. 

Clearly,  anharmonicity  must  be  considered  when  calculating  the 
lattice  contribution  to  the  specific  heat  from  the  phonon  density 
of  states  determined  from  neutron  diffraction  data.  The  specific 
heat  of  the  charge  carriers  is  similar  to  the  predictions  for  a 
Fermi  liquid,  but  differs  markedly  from  the  temperature 
dependence  expected  of  a  gas  of  charged  Bosons.  Also,  the 
experimentally  determined  electronic  density  of  states  is  a 
factor  of  three  larger  than  the  value  predicted  by  band  structure 
calculations.  The  sensitivity  of  many  measured  electronic 
properties  to  the  structural  integrity  of  the  chains  can  be 
attributed  to  the  large  fraction  of  the  density  of  states  at  the 
Fermi  level  contributed  by  the  chain  band. 
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function  for  estimating  the  phonon  density  of  states  from  the 
neutron  density  of  states.  We  thank  V.Z.  Kresin,  J.L.  Cohn,  G. 
Deutscher,  P.B.  Allen,  W.E.  Pickett,  and  R.A.  MacDonald  for  many 
helpful  discussions.  We  credit  R.J.  Soulen,  Jr.  for  stimulating 
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Figure  Captions 

Figure  1.  Relative  enthalpy  of  YBa2Cu307_^  between  t  and  O'C.  The 
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solid  line  through  the  data  represents  a  polynoinial  fit 
to  the  data.  See  discussion  preceding  Eg.  1. 

Figure  2.  The  specific  heat  at  constant  volume  vs.  temperature 

for  YBa2CU307.2,  derived  from  the  measured  enthalpy  data. 
See  Eg.  1.  The  small  circles  are  the  low  tem.perature 
data  measured  by  Junod’^  et  al.  The  lines  represent 
the  lattice  specific  heat  determined  by  Eg.  2  (solid) 
and  Eg.  3  (dotted) . 

Figure  3.  The  electronic  specific  heat  vs.  temperature  for 

YBa2Cu307.2.  In  the  inset,  the  lattice  is  assumed  to  be 
totally  harmonic.  For  the  main  graph,  the  lattice 
contribution  has  been  adjusted  for  anharmonic  effects. 
The  open  and  closed  circles  are  defined  as  in  Fig.  2. 
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